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(N ! Abstract 



By means of modified extended direct algebraic method (MEDA) 
the multiple exact complex solutions of some different kinds of non- 
■ linear partial differential equations are presented and implemented in 

a computer algebraic system. New complex solutions for nonlinear 
equations such as the variant of the RLW equation, the variant of the 
PHI-four equation and the variant Boussinesq equations are obtained 
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^ ■ 1 Introduction 

It is well known that the nonlinear physical phenomena are related to 
■ nonlinear partial differential equations which are involved in many fields 

£NJ . such as physics, chemistry, mechanics, etc. As mathematical models of the 

phenomena, the investigation of exact solutions of the partial differential 
equations will help one to understand these phenomena better. In recent 
years, many powerful methods to construct exact solutions of nonlinear par- 
tial differntial equations have been established and developed. Among these 
are variational iteration method [3l [TJ [2j (6J H], tanh function method 
[2 El, modified extended tanh function method [LQl Q2J Q3J El EE], sine- 
cosine method [161 E] Exp-method[18| , inverse scattering method |19| . Hi- 
rota's bilinear method [20jthe homogeneous balance method [2T]. the Riccati 
expansion method with constant coefficients [23^ [22] . The regularized long 
wave (RLW) equation is an important nonlinear wave equation. Solitary 
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waves are wave packet or pulses, which propagate in nonlinear dispersive 
media. Due to dynamical balance between the non-linear and dispersive 
effects these waves retain a stable wave form. The regularized long wave 
(RLW) equation is an alternative description of non-linear dispersive waves 
to the more usual Kortewege-de vries (KDV) equation [24 1 . Numerical solu- 
tions based on finite difference techniques [251 EE] • Rung-Kutta method [27] 
and Galerkin's method [28] have been given. Alexander and Morris [2H] 
constructed a global trial function mainly from cubic splines. Gardner and 
Gardner [29], using the Galerkin's method and cubic B-spline as element 
shape function to construct an implicit finite element solution. The least 
squares method using linear space-time finite elements used to solve the 
RLW equation [30], Soliman and Raslan [3TJ solved the RLW equation by 
using collocation method using quadratic B-spline at the mid point. Soliman 
and Hussien [32] used the collocation method with septic spline to solve the 
RLW equation, Soliman [2] using the finite difference method with the sim- 
ilarity solution of the partial differential equations to obtain the numerical 
scheme for RLW equation this approach eliminate the difficult associated to 
the boundary, in the end Soliman ]2] using the variational iteration method 
to obtain the exact solutions of the gneralized RLW. The PHI-four equation 
is considered as a particular form of the Klein-Gordon equation that model 
phenomenon in particle physics where kink and anti-kink solitary waves 
interact [33] . Recently, the direct algebraic method and symbolic computa- 
tion have been suggested to obtain the exact complex solutions of nonlinear 
partial differential equations [HTJ [351 [36]. The aim of this paper is to extend 
the modified extended direct algebraic method MEDA method to solve three 
different types of nonlinear differential equations such as the variant of the 
RLW, the variant of the PHI-four and variant Boussinesq equations [12\ \TT\ . 



2 Modified extended direct algebraic method 



To illustrate the basic concepts of the modified extended direct algebraic 
(MEDA) method. We consider a given PDE in two independent variables 
given by 

F(u,u x ,u t ,u xx , ....) = (1) 

We first consider its travelling solutions u(x,t) = u(z), z = i(x + ct) 
or z = i{x — ct), i = \f— 1, then Eq. ([T]) becomes an ordinary differential 
equation 

H{u,iu',-icu',-u",...) =0 (2) 
where v! = j^. In order to seek the solutions of Eq. ([T]), we introduce the 
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following ansatze 

M 

u(z) = a + ^(aj cp 3 ' + bj <p~ j ) (3) 

3=1 

<t>' = b + 4? (4) 

where b is a parameter to be determined, <j) = 'K- 2 )) $ = 757- The 
parameter M can be found by balancing the highest-order derivative term 
with the nonlinear terms [38j . Substituting Q into (J2J) with will yield a 
system of algebraic equations with respect to aj, bj, b and c(where j = 1...M) 
because all the coefficients of have to vanish, we can then determine ao, 
aj, bj, b, and c. Eq. (|3J) has the general solutions: 

(I) If b < 

4> = —\/—b tanh (V—b z), or <f> = —y—b coth (V—b z) 
it depends on initial conditions. 
(II) If b > 

4> = \fb tan (Vb z), or <p = —Vb cot (Vb z) 
it depends on initial conditions. 
(Ill) If 6 = 

^ (5) 

Substituting the results into d3|), then we obtain the exact travelling wave 
solutions of Eq. ([1]). To illustrate the procedure, three examples related to 
variant of the RLW, variant of the PHI-four and variant Boussinesq equa- 
tions are given in the following. 

3 Applications 

3.1 Variant of the regularized long- wave (RLW) equation 



Let us first consider the nonlinear variant of the regularized long wave equa- 
tion which has the form 

ut + a u x - A (u n ) x + (u n ) xxt = (6) 

where a, A, and /? are arbitrary constants. In order to solve Eq. (0) by 
the MEDA method, we use the wave transformation u(x, t) = U (z) with 
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wave complex variable z = i(x — ct), Eq. ([6]) takes the form of an ordinary 
differential equation as 

{a — c)U' — \ (U n )' + c {U n )'" = (7) 

Integrating Eq. ([7]) once with respect to z and setting the constant of 
integration to be zero, we obtain 

(a - c) U - A U n + $ c (U n )" = (8) 

Or equivalently 



(a-c) U- XU n + /3cnU n - L U" + f3cn(n- 1) U n ~ 2 {U'Y = (9) 

Balancing the order of U n with the order of U n ~ l U" in Eq. Q, we find 

M = — ^zj- To get a closed form analytic solution, the parmeter M should 

i 

be an integer. A transformation formula U = V n ~ 1 should be used to 
obtain this analytic solution. So Eq. Q takes the form 

{a-c)(n-l) 2 V 3 -\(n-l) 2 V 2 -(3cn{n-l)VV"+Pcn(2n-l)(V') 2 = (10) 

Balancing the order of V 3 with the order of V V" in Eq. (|10p . gives 
M = 2. So the solution takes the form 

V(z) = a + ai <f>(z) + a 2 <j>{z) 2 + h ^(z)' 1 + b 2 ct> (z)~ 2 (11) 

Inserting Eq. (jlip into Eq. ()10p and making use of Eq. ([3]), using the 
Maple Package, we get a system of algebraic equations, for ao, ai, a2, &i> ^2 
and 6. We solve the obtained system of algebraic equations give the following 
three cases: 

Case (I): consider, A = (— nX— \+2anao) then, a\ = 0, a 2 = 0, b\ = 0, 

, ok A (n 2 — 2 n+1) A z. a A(n 2 -2n+l) .,, , - 

& 2 = 2n/3A . C = 27T^> 6 = ^TjA ' Wlth a bem g an 

arbitrary constant. 



T w \ o.i 2an\(n 2 — 2n + 1) , A ... 

F(x,t) = a (l + cot (-W — — — i x - t 12 

2 y npA 2 n ao 

so, the travelling wave solution is given by 



9 , % / an A (ra 2 — 2n + 1) , 1 A . , . , _^±_ 

u(x, t) = a 1 + cot 2 - W 1 ^ x - t - 1 13 

2 U n p A 2 n ao 
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Case (II): consider, E = ^ "2-2 ^+2^+2/^+02) ^ ^ = 1 a_e ^ 

d = 0, 6i=0, 6 2 = 0, b=\ a ^ ni , c = a2 (w ~ 2 w+1) , with a 2 being an 
arbitrary constant, then 



IEX , 0/ i I E 8 , ao (n 2 — 2n + 1) ... . . 

v <**> = w (1 - tan Wi^b (I " — !))) (14) 

so, the travelling wave solution is given by 



1 EB 2/ i I EB a 2 {n 2 -2n + l) 



«(*,*) = (1 " tan^(- ^ ^(x - - E ' *))))»" 1 (15) 



Case (III) : a = § «i = 0, 61 = 0, b 2 = ^^f^, 

i 2 -2 

16 /3 ri 2 a 2 ' ° £ 



^ = Tfi~R^^~ — ' c = ° 2 ( n 2ra+i) w ^jj fl2 b e j n g an arbitrary constant, then 



Tjr/ s 1 £A ,„ 1 ,,i / EX , a 2 (n 2 -2n + l) . . 



1 ,,i EX , a 2 (n 2 — 2n + 1) , , , . , 



so, the travelling wave solution is given by 



. . . 1 BE 1 2 .i / £ A , a 2 (n 2 -2n + l) , N 



1 9/ i I EX , a 2 (n 2 — 2n + 1) ,,, S ^I. 
-cot 2 -, — 1--^ >-t ) »-i 17 

2 4 \ p n z a 2 E 



3.2 Variant of the PHI-four equation 



A second important example is the Variant of the PHI-four equation [llj . 
which can be written as 

u u - au xx - Xu + 8 u n = n > (18) 

where a, X, and /? are arbitrary constants. In order to solve Eq. (|18p 
by the MEDA method, we use the wave transformation u(x,t) = U(z) with 
wave variable z = i(x — ct) Eq. ([THJ) takes the form of an ordinary differential 
equation 
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c 2 U" - a U" - A U + p U n = (19) 

Or equivalently 

- A U + (3 U n - (c 2 - a) U" = (20) 

Balancing the order of U n with the order of U" in Eq. ()20p , we find M = 

. To obtain a closed form of the analytic solution we use a transformation 
i 

formula U = V™- 1 , that transforms Eq. (|20|) to 



-A(n-l) 2 F 2 +/3(n-l) 2 1/ 3 + (a-c 2 )(r^-l)Fy" + (a-c 2 )(2-n)(y , ) 2 = 

(21) 

Balancing the order of V s with the order of V V" in Eq. (|2ip . we find 
M = 2. So the solution takes the form 

V(z) = a + ai 0(z) + a 2 <K^) 2 + h (^(z)' 1 + b 2 <t>{z)~ 2 (22) 

Substituting Eq. ([25]) into Eq. (|2"Tj) and making use of Eq. , we obtain 
a system of algebraic equations, for ao, ai, 02j &i, and 6. We solve the 
obtained system of algebraic equations give the following three cases: 

n m A(n+1) n n . n , (n+1) (n 2 -2n+l) A 2 

Case (I): a = \ p a\ = 0, a 2 = 0, 6i = 0, 6 2 = 1 ^ (c a_ tt) — , 

6 = A jTxn^ with c being an arbitrary constant, then 



, . A (n + 1) A (n + 1) 2 /« /A (n 2 - 2n + 1) , \ . . 

"''■"'V'V^^V c»-a (I - C "j (23) 

so, the travelling wave solution is given by 

\\(^t) 



, , /A(n+1) A(n + 1) „/» /A(n 2 -2n + l). 



(24) 



/-i /ttN A(n+1) n 2(—an+c 2 n—a+c 2 ) i n , n 

Case (II): a = \ p ai = 0, a 2 = /3( w 2 -2n+i) — ' bl = °' ° 2 = °' 
6 = A 47gz~^) i with c being an arbitrary constant, then 



1 A (n + 1 
2 



1 A 
+ 2- 



-a n + c n — a + c 



tan 2 ( 



/3 

/ /A (n 2 



/3(c 2 



a) 



2n + 1 



a 



(x -ct)) 



(25) 
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so, the travelling wave solution is given by 



1 A (n + 1) 1 A (-a n + c 2 n - a + (?) 
= -p +2 ( C 2 _ ~) 

9 .% I X (n 2 — 2 n + 1) . , \ 

tan ( 2 V c 2 -a (X " Ct)))n " (26) 



Case(III): Consider, L> = ^Xn-Xn^-x-iGab ^ aQ = 1 M^hi) > Ql = , 

° 2 = I A ' &i = 0, 62 = -g ^ TO+ ^ A fe ; c = | -D with 6 being an arbitrary 
constant, then, 



A (n + 1) A (n + 1) 9 . , — - . z „ 
V(x, t) = l 4/j ' + l 8/j ; tan 2 (^ { X --Dt)) + 

^^ C ot\^b( X - l -Dt)) (27) 

so, the travelling wave solution is given by 

. . .A (n + l) A (n + 1) 9/ / — - , i _ 
*) = ( 4/g + V 8/3 ^ tan 2 (^ (z - - D t)) + 

^|±^cot 2 (^(x-^t)))^ (28) 



r< mn 1 A (n+1) n 1 (n+1) A , n , i (n+l) A & 

Case(IV): a = 5 ^ ; , ai = 0, a 2 = 5 v ^ , 61 = 0, 6 2 = g 1 — ^ — , 
c = -7D with 6 being an arbitrary constant then, 



A (n + l) A (n + 1) 9 . , — - , i „ .. 
V(x, t) = l 4/j } + l 8/j ; tan 2 (^ (x + -D t)) + 

^^ CO t\^b{ X+ l -Dt)) (29) 

so, the travelling wave solution is given by 

. . .A (n + l) A (n + 1) 9/ / — r . i _ 
t} = ( 4/ 3 + 8/3 tan 2 (^ (z + - D t)) + 

^|±^cot 2 (^6(x + ^t)))^ (30) 
All the solutions of the equations are new. 
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3.3 The variant Boussinesq equations 



Finally, we consider a very important example as an illustration of the mod- 
ified extended direct algebraic method for solving the variant Boussinesq 
equations |12j . we will consider the following system of equations 

u t + v x + u u x = (31) 
v t + (uv) x + u xxx = (32) 

To solve the system of Eqs (|3T|) . ([32]) by means of the modified 
extended direct algebraic method, we use the wave transformation u(x, t) = 
U (z) and v(x, t) = V{z) with complex wave variable z = i(x+Xt). Therefore, 
system (f3Tj) , (f32|) is reduced to the ordinary differential equations in the form 

XU' + V' + ^ (U'f = (33) 
XV' + {UV)' - U'" = (34) 

Integrating both equations once leads to: 



Ci - XU ~\u 2 = V (35) 
XV + U V - U" = C 2 (36) 

where C\ and C2 are integrating constants, so as to we find the special 
forms of the exact solutions, for simplicity purpose, we take C\ =Ci = 
Substituting ([35]) into l[36|) gives 



U" + ^U 3 + ^XU 2 + X 2 U = 0, (37) 

By balancing U" with U ?J in eq. (|37|) . we find M = 1. So the solutions 
take the form 

U(z) = a + a 1( f)(z) + b l( f)(z)- 1 (38) 



V(z) = -X [a + a x <t>{z) + h^z)' 1 ] - i [a + a X (j>{z) + h^z)' 1 ] 2 (39) 

Inserting Eqs. ([38]) . (|4T|) into Eq. ([36]) . making use of Eq. (jlj), and by 
using Maple Package, we get a system of algebraic equations, for ao, ai, &i, A, 
and b. We solve the obtained algebraic system of equations by Mable Pack- 
age and select four cases of solutions as: 
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Case (I): a\ = 2i, b\ = 0, b = a 2 ,, A = — ao, with ao being an 
arbitrary constant, the complex wave solutions are : 

u(x, t) = ao(l + ztan(— (x — aot))) (40) 



v(x,t) = ao(ao(l + *tan(y(x-a t))))-^(ao(l + «tan(y(a:-aot)))) 2 (41) 



Case (II): a\ = — 2i, b\ = 0, b = a 2 ,, A = — ao, with ao being an 
arbitrary constant, the complex wave solutions are : 

u(x,t) = ao(l — ztan(— (x — aot))) (42) 



v(x,t) = ao(ao(l-jtan(y(x-a t))))-^(a (l-itan(^(a:-aot)))) 2 (43) 



Case (III): a\ = 2i, b± = a\i, b = a 2 ,, A = — ao with ao being an 
arbitrary constant, the complex wave solutions are 

u(x, t) = a + V ^ a ° ( tan ( ^ Q ° (x - a t)) + cot ( ^ Q ° (x - a t))) (44) 



v(x,t) = a (a + V ^ Q ° (tan ( [ x - a t)) + cot ( (x - a t)))) 

-^(°o(ao + (tan ( ^°° (x - a t)) + cot ( ^ Q ° ( x _ ao t)))) 2 (45) 

Case (IV):: a± = — 2 i, b\ = \ a\ i, b = | a 2 ,, A = — ao with ao being an 
arbitrary constant, the complex wave solution are 



u(x, t) = ao — ° (tan ( (x — ao i)) + cot ( - (x — ao £))) (46) 
w(x,t) = a (a - V ^ a ° (tan ( V ^ Q ° ( x - a i)) + cot ( (x - a i)))) 

-^ao - 4^( tan - °° t)] + cot - a ° t))))2(47) 

All the solutions of the variant Boussinesq equations are new. 
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4 Conclusions 



In this paper, the MEDA method has been successfully applied to find the 
solution for three nonlinear partial differential equations such as the the 
variant of the RLW equation, the variant of the PHI-four equation, and 
the variant Boussinesq equations . The modified extended direct algebraic 
method is used to find a new complex travelling wave solutions. The re- 
sults show that the modified extended direct algebraic method is a powerful 
mathematical tool to solve the the variant of the RLW equation, the variant 
of the PHI-four equation, and the variant Boussinesq equations it is also a 
promising method to solve other nonlinear equations. 
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